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The decay Kl —t invisible has never been experimentally tested. In the Standard Model (SM) 
its branching ratio for the decay into two neutrinos is helicity suppressed and predicted to be 
Bc{Kl —» vv) < 1CP 10 . We consider several natural extensions of the SM, such as two-Higgs-doublet 
(2HDM), 2HDM and light scalar, and mirror dark matter models, those main feature is that they 
allow to avoid the helicity suppression factor and lead to an enhanced Bt(Kl —» invisible). For the 
decay Kl -A vv the smallness of the neutrino mass in the considered 2HDM model is explained by 
the smallness of the second Higgs doublet vacuum expectation value. The small nonzero value of 
the second Higgs isodoublet can arise as a consequence of nonzero quark condensate. We show that 
taking into account the most stringent constraints from the K —>• 7r + invisible decay, this process 
could be in the region of Bt(Kl —> invisible ) ~ 10 -8 — 10 -6 , which is experimentally accessible. 
In some scenarios the Kl -A invisible decay could still be allowed while the K —¥ n + invisible 
decay is forbidden. The results obtained show that the Kl -A invisible decay is a clean probe of 
new physics scales well above 100 TeV, that is complementary to rare K —> n + invisible decay, and 
provide a strong motivation for its sensitive search in a near future experiment. 


PACS numbers: 14.80.-j, 12.60.-i, 13.20.Cz, 13.35.Hb 


I. INTRODUCTION 

In the Standard Model (SM) the branching ratios of 
the I\ + -A 7r + + invisible and Kl —> 7r° + invisible 
decays are predicted to be jjJ 

Br(K L -A t t°vv) = (2.6 ± 0.4) x 10" 11 , (1) 

Br(K + 7 t + vv) = (8.5 ± 0.7) • 1CT 11 , (2) 

with the invisible final state represented by neutrino 
pairs. A strong comparison between experiment and the¬ 
ory is possible due to the accuracy of both the measure¬ 
ments and the SM calculations of these observables. A 
discrepancy would signal the presence of physics beyond 
the Standard Model (BSM) making the precision mea¬ 
surements of these decays an effective probe to search for 
it, see e.g. dj-Q. 

The branching ratio of the Kl —> invisible decay in 
the SM is predicted to be very small compared to those 
of Eq.JI]) and Eq.© for v masses laying in the sub-eV 
region favored by observations of v oscillations jjj}. In¬ 
deed, the Kl has zero spin, and it cannot decay into two 
massless neutrinos, as it contradicts to momentum and 
angular momentum conservation simultaneously. For the 
case of massive vs their spins in the Kl rest frame must 
be opposite and, therefore, one of them is forced to have 
the ’’wrong” helicity. This results in the Kl —> vv 
decay rate being proportional to the v mass squared 

T (K l —> vv) oc < 10 -17 assuming m v < 1 

eV. However, if one take the direct experimental upper 
limit on the v T mass m„ T < 18.2 MeV [7|, the predicted 
branching ratio, calculated at the quantum loop level is 

a 

Br(K L -A W) ~ 1CT 10 (3) 


Therefore, an observed Bt(Kl -a invisible) 10 -10 
would unambiguously signal the presence of BSM 
physics. 

The decay Kl — > invisible has never been experimen¬ 
tally tested. Since long ago it was recognized that this 
decay ’’would be interesting to explore, but its detec¬ 
tion looks essentially impossible. New ingenious exper¬ 
imental ideas are required” [8j. Recently, an approach 
for performing such kind of experiments by using the 
K + n —> K°p (or K~p -A K n ) charge-exchange reac¬ 
tion as a source of well tagged A'°’s has been reported 
[S]. At the same time the first experimental bound 
Bv(Kl — > invisible) < 6.3 • 10~ 4 has been set from ex¬ 
isting experimental data. It has been shown, that com¬ 
pared to this limit, the expected sensitivity of the pro¬ 
posed search is at least two orders of magnitude higher - 
Bt(Kl — > invisible ) < 10” 6 per ~ 10 12 incident kaons. 
It could be further improved by utilizing a more care¬ 
ful design of the experiment, thus making the region 
Bv(Kl —> invisible) ~ 10 -8 — 10 -6 , or even below, ex¬ 
perimentally accessible @. 

Being motivated by these considerations, we discuss 
in this work several natural extensions of the SM and 
show that taking into account the most stringent con¬ 
straints from the measured K + —> n + + invisible de¬ 
cay rate, the decay Kl —> invisible could occur at the 
level Bt(Kl — > invisible) ~ 10 -8 — 10~ 6 . The main 
feature of the considered models, that leads to an en¬ 
hanced branching ratio for Kl —> invisible, compared 
to K + —> 7r + + invisible, is that they allow to avoid 

the helicity suppression factor ^ in the SM, while 

profiting from its larger phase-space due to the decay 
into two light weakly interacting particles. In addition, 
there might be the case when Kl —> invisible could still 
be kinematically allowed, while K + —> 7r + + invisible is 
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forbidden. Additional motivation to search for the AY 
(and Ks) invisible decay is related to precision tests of 
the K° — TG° system by using the Bell-Steinberger uni- 
tarity relation [9(. This relation connects CP and CPT 
violation in the mass matrix to CP and CPT violation 
in all decay channels of neutral kaons and is a powerful 
tool for testing CPT invariance with neutral kaons S3- 
The question of how much the invisible decays of K$ 
or AY can influence the precision of the Bell-Steinberger 
analysis still remains open |lTj | . All this makes the fu¬ 
ture searches for this decay mode very interesting and 
complementary to the study of the K —> 7r + invisible 
decays. 


II. K l -a vv DECAY IN MODEL WITH 
ADDITIONAL SCALAR DOUBLET 

Consider now the Kl, —> vv decay in the two-Higgs- 
doublet model (2HDM) with an additional heavy Higgs 
doublet H 2 • This type of 2HDM models can introduce 
flavor-changing neutral currents, provide explanations of 
the orig in of Dark Matter and CP violation , see e.g. 
Ref. [12j . The interaction of the heavy isodoublet field lf 2 
with quarks, leptons and the standard Higgs isodoublet 
H leading to the Kl -a vv decay has the form 

Pint b , 2TprP^2^ / Tji T h2dLS R Q\Lhi2^R (4) 
-\-5itilljl 2 H+ H 2 + h.c. — A' 1 ^j 2 H 2 H 2 > 

where L T = (v tl ,t l ), Qil = (u L ,d L ), H 2 = 

H 2 = ((H°)*,-(H+)*) and /i 2r , h 2 d L s R are Yukawa cou¬ 
pling constants. Note that in general the second Higgs 
isodoublet H 2 will have nonzero Yukawa interactions 
with other quark and lepton fields but since we are in¬ 
terested mainly in the Kl -a vv decay we have writ¬ 
ten explicitly only the Yukawa interactions important for 
us. In considered model the neutrinos acquire nonzero 
Dirac masses m Vr = /i 2t < H 2 > due to nonzero vac¬ 
uum expectation value of the second Higgs isodoublet H 2 

< H 2 >= ' "," 2 < H > (< H >= 174 GeV) and the 

H 2 

smallnes of the Dirac neutrino masses is a consequence 
of the < H 2 > smallnes. The smallnes of < H 2 > is due 
to the assumed large value of Mjl 2 or(and) small value 
of Smj 1H2 [35]. For instance, for m Vr =0.1 eV, /i 2t = 0.1 

and Mh 2 = 10 5 GeV we find that = 0.6 • 10~ n 

h 2 

and 8m 2 HH = 0.06 GeV 2 . It is interesting to note that 
for 5m 2 HH2 = 0 the < H 2 >= 0 at classical level but 
the spontaneous symmetry breaking of SUl{3 ) <8> SUr(3) 
chiral symmetry in QCD leads to nonzero vacuum expec¬ 
tation values for the Higgs fields [l4[. Really, for monzero 
Yukawa interaction L H2 Q ld = h 2 d L d R QiLH 2 dR + h.c . 
due to nonzero vacuum expectation value of quark con¬ 
densate < dd >= — (/n- = 93 MeV ) the field 

< H 2 > acquires monzero vacuum expectation value < 
H 2 >= 2 h 2 d^d R Ml ■ Numerically for h 2T = h 2dL d R = 1 


and ?riy T = 0.1 eV we find that Mh 2 ~ O(10 4 ) GeV. So 
in this model with 8m 2 HH2 = 0 the vacuum expectation 
value < H 2 >= 0 at tree level but the nonzero quark con¬ 
densate leads to the appearance of small vacuum expec¬ 
tation value < H 2 >/ 0 for the second Higgs isodoublet 
that explains the smallnes of the neutrino masses. 

For the case of nonzero neutrino Majorana mass ni i/TR 
we assume that the mass m,, is small so the decay 
Kl — > v T v T is kinematically allowed. Again, as in the 
previous case we assume that the Dirac neutrino mass 
arises due to nonzero < H 2 > vacuum expectation value 

m 2 

and the smallness of the see saw = D,/T neutrino 

TR m ^ R 

mass is again explained due to the smallness of < H 2 >■ 
The Lagrangian (4) contains AS = 1 neutral flavour 
changing terms but for heavy doublet i? 2 it is not dan¬ 
gerous. The effective four fermion Lagrangian describing 
the decay Kl —> v T v T has the form 

Le ff = jp dLS R v TL v TR + h.c ., (5) 

X 


where 


1 

Ml 


T'2dLSR'T'2T 

Mr 2 


( 6 ) 


As it has been mentioned before we assume the exis¬ 
tence of small Dirac or Majorana neutrino mass v T . The 
decay rate of the invisible decay Kl —> v T v T is deter¬ 
mined by formula 


F (Kl — t v L tVrt,vr t vlt) 


I67 tM^ 


Fk 


2 (m d + m s ) 


fK{m 2 JM 2 KL ), 


(7) 


where K{x) = (1 — 4a:) 1 / 2 for Dirac neutrino with a mass 
m, l/T and K(x) = (1—a:) 2 for Majorana neutrino v TR with 
a mass w„ tr . Here Fk ss 160 MeV is kaon decay con¬ 
stant and m s ,mj are the masses of s- and d-quarks[36j. 
For BrjKr, —>■ v T v T ) = 10 -6 we can test the value of Mx 

up to m 

M x < 0.6 • 10 5 GeV (8) 


for small Dirac or Majorana neutrino mass m Vr <C Mk l ■ 
It should be noted that the existence of AS = 1 neutral 
flavour changing interaction (5) leads to additional con¬ 
tribution to rare decays Kl —> n a vv and K + —> ir + vv. 
The current experimental values are [n|, [l6j 

Hr (AY -A 7 t°vv) < 2.6 x 10 -8 , (9) 


Br{I<+ -A 7 t+vv) = (17.3±io.i) ' 10_11 . (10) 

with the SM predictions of m and ([2|>. respectively. The 
measured value (10) for the Br(K + -a tt + vv) allows to 
set more stringent constraints. Therefore, we restrict our¬ 
selves to the calculation of the BSM contribution only to 
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this decay channel by using the effective Lagrangian (5). 
This leads to the following formula for the differential 
K + —> n + vv decay width: 


dr 


BSM 


( K + —> TT + w) 1 


dq 2 


(g 2 - 


V Tj r/ 


(2tt) 3 32 M 3 k+ 


q 2 M* 

‘ \J+ M 2 + - g 2 ) 2 - 4M| + M 2 4 
/ 0 (g 2 )(M£ + -m 2 + ) 


2 {-m d + m s ) 


( 11 ) 


The form factor /o(g 2 ) is determined in standard way 


as 


< n\d^s\K >= f + (q 2 ){P K +P„y+f_(q 2 )(P K -P„y = 

( 12 ) 


f + (q 2 )[(P K + P„r~ 


Ml - Ml 


q 


q^ + foiq 2 ) 


M 2 K - Ml 


q 


-q 


where g M = (Pk — P-kY and m 2 R < ( M K + — M^+) 2 . The 
form factors /+ and /o are related to the exchange of 1“ 
and 0+, respectively. The following relation holds: 

/+(0) = /o(0) , /o(g 2 ) = f+(q 2 ) + u o q _ 1 / 2 /_(g 2 ). 

(13) 

In our calculations we use standard linear parametriza- 
tion for the form factor / 0 (g 2 ), namely 

/o(g 2 )=/o(0)(l + A o -|J-). (14) 

Numerically we take /o(0) = 0.96(11} and Ao = —0.06 

Gi¬ 
lt is convenient to represent the result in terms of the 

ratio fi~ l = because the ratio /? does not 

depend on unknown value of Mx- Also p does not de¬ 
pend on the values of quark masses md, m s . For the case 
of massless neutrino we find that 


/3 ~ 2 • 10 -3 . (15) 

Note that the smallness of the p is mainly due to the 
3-body phase space smallnes in comparison with 2-body 
phase space. From the difference between the theoretical 
and experimental values (2) and (10), respectively, by 
summing up errors of (10) in quadrature we find that 
the BSM contribution to the Br(K + — >■ 7 t + vv) is less 
than 


Br BSM (K+ -A n+vv) < 2.1 • 1CT 10 . (16) 

From the limit (16) and the estimate (15) we find that 
for massless neutrinos 

Br{K L vv) < 1CT 7 (17) 

The estimates (15, 17) are valid for small m VR -C M n + 
Majorana msss of righthanded neutrino. For higher m UR 


values the limit (17) is more weak and for the case Mk l > 
m„ B > M k + — M n + when the decay K + —> n + v TL v TR is 
kinematically prohibited, but the decay K R —> vv is still 
allowed, the restriction from K + —> n + vv decay does not 
work. 

The measured (K R — K R ) mass difference strongly re¬ 
stricts [2(| the effective AS = 2 interaction 

Lsdsd. = to— s~ R dLS~ndL + h.c. . (18) 

A sdsd 

Namely l2tl 


Asdsd A 1.8-10 7 GeV. (19) 

For the model (2) with the the additional Higgs doublet 
Hi = {H 2 , 772,1 + *# 2 , 2 ) we find that 


A 2 

sdsd 


— 1^2 d.L& 


\h 2 d L s R \ 2 ^ m HH 2 


Mlo 

ri 2,\ 


M h° 

n 2,2 


Mk 


M'k 


( 20 ) 

Using the bound (19) we can restrict the parameter 
8m 2 H H 2 ■ For instance, for Mh 2 = 10 5 GeV, h 2 d LSR = 1 
we find 8m 2 HH2 < 0.3 • 10 6 GeV 2 that is much more weak 
than the estimate of 8m 2 HH2 coming from the neutrino 
mass. 

In general case we can have additional flavour chang¬ 
ing Yukawa interaction h 2sL d R Q 2 LH 2 d R + h.c (Q 2 l = 
( c L ,s L ) in the Lagrangian (4) that leads to the tree level 
flavour changing AS = 2 effective interaction L e // = 


2d ^ 2 J LdR ( dLS R d R SL + h.c.) We can simultaneously 

avoid the AS = 2 bound A As =2 = {h 2dL s R h* 2sLdR Y 1/2 ■ 
Mh 2 > 1.8 • 10 7 GeV and obtain phenomenologically 
interesting values for Br(K R —> vv) for small quark 
Yukawa coupling constants h 2 d LSR , h 2sL d R , relatively 
light second Higgs doublet and not small lepton Yukawa 
coupling constant h 2r . For instance, for h 2 d LSR = 
h 2sL d R = (1/300) 2 , h 2r = 1 and Mh 2 = 300 GeV we 
find that A As =2 = 2.7 • 10 7 GeV and Br(K R —> vv) = 
0.4 • 10 -6 . The existence of relatively light with a mass 
Mh 2 = 300 GeV second Higgs doublet does not contra¬ 
dict the LHC data. The best way to look for the second 
Higgs isodoublet at the LHC is the use of the reaction 
pp —> Z* / gamma* —> Hfi —>• t + t~vv. So the sig¬ 
nature is two r leptons plus nonzero EP SS in final state 
that coincides with the signature used for the search for 
direct production of stau leptons at the LHC. 


III. K l -j- </></> DECAY IN MODEL WITH 
ADDITIONAL SCALAR DOUBLET AND 
SCALAR SINGLET 0 


Consider now the K R —> (fxf> decay in the extension 
of the SM with heavy Higgs doublet H 2 and light neu¬ 
tral scalar singlet field 4>. The Yukawa interaction of 
the heavy isodoublet H 2 with quarks and the interaction 
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of the <j> field with Higgs isodoublets H 2 and H{ Higgs 
isodoublet of the SM) has the form 

Lr = h 2dL s r QilsrH 2 + A (H+H)4> 2 + 6mjj H2 H+H{ 21) 

+h.c. — Mfi 2 H 2 H 2 , 

where Q 1L = ( u L ,d L ), H 2 = and h 2dLSR , 

A are Yukawa and Higgs couplings. After electroweak 
SUr(2)®U(1) symmetry breaking trilinear term describ¬ 
ing transition H 2 —> (jxj) 


LH 2 if><f> — A < H > H 2 (f>~ + h.c. 
arises. The effective Lagrangian 

1 - 2 

L e ff = -j^dLSRct) + h.c ., 

1 _ h 2dLSR X <c H 7 

~Mx ~~ Mf a 


( 22 ) 


(23) 


(24) 


describes invisible decay Kr —> (jxj). Here we assume that 
the mass of (j> is less than Mr l /2. The decay rate of the 
invisible decay Kr —> (jxj) is determined by formula 

lVf 3 rp 

r ^ - +*> = ’ 

(25) 

where K(x) = (1 — 4a;) 1 / 2 . For Br{Kr —»• (jxj)) = 10 6 
and <C Mr l we can test the value of Mx up to 


M x < 10 10 GeV. 


(26) 


For A = 1 and h 2 d LSR = 1 the mass of the second Higgs 
isodoublet can be tested up to Mr 2 < 10 6 GeV. 

The bound (16) allows to restrict the Kr —> (jxj) de¬ 
cay in full analogy with previous model. Namely, in the 
model with the effective Lagrangian (22) the Kr —> (jxj) 
decay width is determined by the expression 


dT BSM ( R + _+ 
dq 2 


1 


1 


(27r) 3 32 M\ l M| 
4 m 2 , 


! - ^±m~L 

'\J[(Mk+ + M 2 + - 9 2 ) 2 - 4M| + M 2 + ](1 - -/) 


\fo{q 2 ){M%+ -m 2 + ) 


2 (~m d + m s ) 


(27) 


It is convenient to use the ratio /? 1 = 
because the ratio /3 does not depend on unknown value 
of Mx and on the values of quark masses m d , m s . For 
the case m# <C M n + we find that 


j3 w 10 


-2 


(28) 


As in the previous model the smallness of the /3 is mainly 
due to the 3-body phase space smallnes in comparison 
with 2-body phase space. 


From the (16) and (28) we find 

Br(K r ->• <j)(j)) < 2 • 1CT 8 . 


(29) 


For not very light ^-particle the limit on Br{Kr —>• <jxj>) 
will be not so stringent as the bound (29), moreover, for <j> 
particle mass Mr l /2 > m$ > (M K +—M 7T +)/2 the decay 
K + —>• tt + dD is kinematically prohibited while the decay 
Kr — > (jxj) is allowed. Therefore the bound (29) derived 
from the decay width of the I\ + —> Tt + vv decay does not 
work for Kr —>• (jxj) decay mode. Note, that such sub-Gev 
scalar (j> could be a good dark matter candidate [ 2 l| . As 
in the previous model the bound from the Kr — Kg mass 
difference leads to the bound on the unknown parameter 
Sm‘j IH 2 the level 5m 2 HH2 < 30 GeV 2 for = 

10 4 GeV. 


IV. K l invisible DECAY IN MODEL WITH 
MIRROR WORLD 

Finally, we discuss the Kr oscillations into a hidden 
sector, which would manifest themselves through the 
Kr —>• invisible decay. As an example of such hidden 
sector we consider the one of the mirror matter models. 
The idea that along with the ordinary matter may exist 
its exact mirror copy, introduced for the parity conserva¬ 
tion, is not new |22|. Accordingly, each ordinary particle 
of the SM has a corresponding mirror partner of exactly 
the same mass as the ordinary one. The mirror fields are 
all singlets under the SM SU C ( 3) ® SUr{2) ® 17(1) gauge 
group. Mirror matter is dark in terms of the SM inter¬ 
actions, and could be a good candidate for dark matter, 
see, e.g., Refs.(23|, and recent 24[. In addition to gravity, 
the interaction between our and this type of dark mat¬ 
ter could be transmitted by some gauge singlet particles 
interacting with both sectors. Any neutral, elementary 
or composite particle, in principle, can have mixing with 
its mirror duplicate. This results in several interesting 
phenomena, such, e.g. as Higgs 25|, positronium j26j . 
muonium [2?[ . or neutron [28j oscillations into their hid¬ 
den partner, which have been or planned to be experi¬ 
mentally tested I29l - f32i] . 

In particular, the neutral A^-meson can mix with it 
mirror (to) analog Kf j m due to effective four-fermion in¬ 
teraction 

Lint = “ 75)ss m 7/i( 1 - 75 )d m ] (30) 

m 

The interaction (30) leads to conversion of ordinary Kr- 
meson to mirror AT,-meson. The decays of mirror Kr- 
meson are invisible in our world that leads to invisible 
Kr decay with the branching ratio 

S 2 

Br(Kr —> invisible) = —^— T - w > (31) 


where 


<5 = 


Mi 


K l 


2( ( 5 2 +^ 2 ot (A i )) , 


^ L ,m\Lint\-^ L ^ • 


(32) 
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For the interaction (27) in the vacuum insertion approx¬ 
imation we find that 


f 2 k m Kl 

Ml 


(33) 


Numerically, for Bv(Kl —>• invisible) = 10 6 we can 
probe the value of M m up to 

M m < 8.4 • 10 s GeV . (34) 

In our estimates we used nonrenormalizable effective 
four-fermion interaction (30). It is possible to obtain 
the effective interaction (27) from renormalizable mirror 
world model with the Higgs doublet extension of the SM 
model (see previous discussions) and with the additional 
interaction term between our and mirror world 


L m = Xm(H+H 2 )(H m Hl 2 ) + h.c. . (35) 

After electroweak symmetry breaking in our and mirror 
worlds (< H >=< H m >rs 174 GeV) we find an effective 
four-fermion interaction 

Leff = -r72^LSRSR^md Lyrn + h.C. , (36) 

lV±m. 


where 


1 = h ld LSR X m \<H> | 2 
Ml M| 2 ' Ml 2 


10~ 10 would unambiguously signal the presence of BSM 
physics. We consider the Kr —> invisible decay in sev¬ 
eral natural extensions of the SM, such as the 2HDM, 
2HDM and light neutral scalar field </>, and mirror dark 
matter model. Using constraints from the experimen¬ 
tal value for the Br(K + — > 7 t + vv) we find that the 
Kl —>• invisible decay branching ratio could be in the re¬ 
gion Bt(Kr —>• invisible ) ~ 10~ 8 — 10~ 6 , which is exper¬ 
imentally accessible allowing to test new-physics scales 
well above 100 TeV. In some scenarios these bounds 
can be avoided, as in the model with the massive right- 
handed neutrino and scalar </>-particle. This makes the 
Kl —» invisible decay a powerful clean probe of new 
physics, that is complementary to other rare K decay 
channels. Additionally, in the case of observation the 
Kl —> invisible decay could influence the precision of 
the Bell-Steinberger analysis of the K° — K system. 
The results obtained provide a strong motivation for a 
sensitive search for this process in a near future K de¬ 
cay experiment proposed in [$]. It should be noted that 
in full analogy with the case of Kl invisible decay we 
can expect the existence of invisible decays of Bd and B s 
mesons, see e.g. (H, (34J, with the branchings similar to 
those discussed above. 


V. CONCLUSION 


In conclusion, the observation of the Kl invisible 
decay with the branching ratio Bv{Kl —> invisible) 


This work was supported by the Grant RFBR N 13- 
02-00363. 
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